Three-body force effect on nucleon momentum distributions in asymmetric nuclear 
matter within the framework of the extended BHF approach 
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We have investigated the three-body force (TBF) effect on the neutron and proton momentum 
distributions in asymmetric nuclear matter within the framework of the extended Brueckner-Hartree- 
Fock approach by adopting the Al/lS two-body interaction plus a microscopic TBF. In asymmetric 
nuclear matter, it is shown that the neutron and proton momentum distributions become different 
from their common distribution in symmetric nuclear matter. The predicted depletion of the proton 
hole states increases while the neutron one decreases as a function of isospin-asymmetry. The TBF 
effect on the neutron and proton momentum distributions turns out to be negligibly weak at low 
densities around and below the normal nuclear density. The TBF effect is found to become sizable 
only at high densities well above the saturation density, and inclusion of the TBF leads to an overall 
enhancement of the depletion of the neutron and proton Fermi seas. 
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^ ■ I. INTRODUCTION 

To determine reliably the properties of isospin asymmetric nuclear matter is a challenge in nuclear physics and 
nuclear astrophysics [lHj|- The nucleon momentum distribution in nuclear matter is one of the most important 
properties of nuclear matter. Many-body correlations induced by nucleon-nucleon (NN) interactions among nucleons 
' q1 play a significant role in a nuclear many-body system, which make the system much more complicated and have 
^ more plentiful properties than a non-interacting Fermi system. For example, the effect of the short-range correlations 
Ch may lead to the depletion of the nucleon momentum distribution below the Fermi momentum and the population 
above the Fermi momentum in nuclear matter [sj . The nucleon momentum distribution is of a great physical interest 

■ since it is closely related to the nature of the underlying NN interaction. The depletion of the Fermi sea is expected 
J> , to be closely related to the hard core and the tensor component of the NN interaction It plays an important 

00 ' role in testing the validity of the physical picture of independent particle motion in the mean field theory or the 

^""^ ' standard shell model and serves as a measure of the strength of the dynamical NN correlations induced by the NN 

• interaction in a nuclear many-body system [t], Q . The study of the nucleon momentum distribution in nuclear matter 

. \ may provide desirable information on the depletion of the deeply bound states inside finite nuclei and is expected 

■ to be important for understanding the structure of finite nuclei. Experimentally, the effects of NN correlations and 
' the nuclear depletion of the Fermi sea can be investigated by the (e,e'p), (e, e'NN), and proton induced knock-out 
, reactions (QhIIi] . The related measurements have been reported continually |T2| - [2^ and definite evidence of the short- 

. . ■ range NN correlations has been observed in these experiments. The analysis of the (e, e'p) reactions on •^°^Pb at 
^ NIKHEF has indicated that a depletion of 15% — 20% for the deeply bound proton states is required for describing the 
measured coincidence cross sections [ISj] . Recent experiments on the two- nucleon knock-out reactions have shown that 
nucleons can form short-range correlated pairs with large relative momentum and small center-of-mass momentum. 
d ■ A strong enhancement of the neutron-proton {np) short-range correlations over the proton-proton (pp) correlations 
has been observed at JLab [l^ due to the dominate role played by the short-range tensor components of the NN 
interactions in generating the NN correlations [23j , and may have significant implications for neutron star physics [24| . 

The short-range correlations and the nucleon momentum distribution in nuclear matter have been investigated 
extensively by using various theoretical ap proaches, such as the extended Brueckner-Hartree-Fock (BHF) method 
[25l - [3l| . the Green function theory [33 - l36l |. the in- medium T-matrix method [13, HI], the variational Monte Carlo 
approach [11], and the correlated basis function approach [s^ The predicted de plet ion of the Fermi sea by 

adopting different theoretical approaches has been shown to be slightly larger than 15% [26l l29ll39l, l40j . In Ref. [29| . 
the authors have calculated the nucleon momentum distribution and quasiparticie strength in symmetric nuclear 
matter in the framework of the Brueckner-Bethe-Goldstone theory by including high-order contributions in the hole- 
line expansion of the mass operator, and a good agreement between the calculated quasiparticie strength and the 
experimental data [l^ has been shown. In Ref. |3l| . the neutron and proton occupation probabilities, averaged 
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below their respective Fermi seas in neutron-rich matter at the saturation density, have been predicted. In more 
recent papers (sBI . |36| , the nucleon momentum distributions in neutron matter and asymmetric nuclear matter have 
been investigated within the framework of the Green function method, and the isospin-asymmetry dependence of the 
depletions of the neutron and proton Fermi seas has been clarified. It has been shown fsij, [s^ [s^ that increasing the 
isospin asymmetry leads to an increasingly larger depletion of the proton hole-states than that of the neutron hole 
states in asymmetric nuclear matter. In the present paper, we shall extend the previous investigation of Ref. j29l | to 
asymmetric nuclear matter and investigate the isospin-asymmetry dependence of the neutron and proton momentum 
distributions within the extended Brueckner-Hartree-Fock (EBHF) method. Particulary, we will concentrate on the 
TBF effect on the momentum distributions and their isospin dependence in asymmetric nuclear matter especially at 
suprasaturation densities. 

The present paper is organized as follows. In the next section, we give a brief review of the adopted theoretical 
approaches including the EBHF theory and the TBF model. In Sec. Ill, the calculated results will be reported and 
discussed. Finally, a summary is given in Sec. IV. 



II. THEORETICAL APPROACHES 

The present calculations are based on the extended BHF approach for asymmetric nuclear matter [4l| . The extension 
of the BHF scheme to include microscopic three-body forces can be found in Refs. (42| - |4^ . Here we simply give a 
brief review for completeness. The starting point of the BH F ap proach is the reaction G-matrix, which satisfies the 
following isospin dependent Bethe-Goldstone (BG) equation p5| : 

n( a \ , 1/ \kik2)Q{ki,k2){kik2\ . . , , 

Z_. ^ _ ,i^k,) - e{k2) 

where ki = {ki, Ui, Ti) denotes the momentum, and the z-component of spin and isospin of a nucleon, respectively. Vmn 
is the realistic NN interaction and w is the starting energy. The asymmetry parameter (3 is defined as /? = (/9„ — Pp)/ p, 
where p, pn and pp denote the total nucleon, neutron, and proton number densities, respectively. The Pauli operator 
is defined as Q(fci,fc2) = [1 — no(fci)][l — fto(^2)]i a-nd it prevents two nucleons in intermediate sates from being 
scattered into their respective Fermi seas (Pauli blocking effect). Here by no(fc) we denote the Fermi distribution 
function which is a step function at zero temperature, i.e., rio(fc) = 9{kp — k). The single-particle (s.p.) energy e(fc) 
is given by: e{k) = h'^k'^ / (2m) + U{k), where the auxiliary s.p. potential UsBFik) controls the convergent rate of the 
hole-line expansion j45j]. In the present calculation, we adopt the continuous choice for the auxiliary potential since 
it provides a much faster convergence of the hole- line expansion up to high densities than the gap choice [i^ . Under 
the continuous choice, the s.p. potential describes physically at the lowest BHF level the nuclear mean field felt by a 
nucleon in nuclear medium (47i] . and is calculated as follows: 

U{k)^Re {kk'\G[p,e{k)+e{k')]\kk')A , (2) 

k'Kkp 

where the subscript A denotes anti-symmetrization of the matrix elements. 

For the realistic NN interaction Vattv, we adopt the Argonne Vis (^V^18) two-body interaction \a§\ plus a micro- 
scopic TBF [4^ constructed by using the meson-exchange current approach In the TBF model adopted here, 
the most important mesons, i.e., tt, p, a, and w have been considered. The parameters of the TBF model, i.e., the 
coupling constants and the form factors, have been self-consistently determined to reproduce the A^IS two-body 
force using the one-boson-exchange potential model and their values can be found in Ref. j43| . In our calculation, the 
TBF contribution has been included by redu cing the TBF to an equivalently effective two-body interaction according 
to the standard scheme as described in Ref. 43. In r-space, the equivalent two-body force V^^ reads 



xW:i{r;r^r^\rAr^)cj>,,{r^){l - r;(ri3))(l - r;(r23)). (3) 

In order to calculate the nucleon momentum distribution in nuclear matter with the EBHF approach, we follow the 
scheme given in Refs. 0, [2^ and extend the scheme to asymmetric nuclear matter. Within the framework of the 
Brucckner-Bethe-Goldstone theory, the mass operator can be expanded in a perturbation series according to the 
number of hole lines, i.e., 

Ar(fc,w) = M{{k,u}) + M^{k,u}) + M^{k,Lo) + ■■■ (4) 
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where t denotes neutron or proton (hereafter we will write out explicitly the isospin index r). The mass operator 
is a complex quantity and the real part of its on-shell value can be identified with the potential energy felt by a 
neutron or a proton in asymmetric nuclear matter. In the expansion of the mass operator, the first-order contribution 
M^{k,uj) corresponds to the standard BHF s.p. potential and the on-shell value of its real part coincides with the 
auxiliary potential under the continuous choice given by Eq. ([2]). The higher-order terms stem from the density 
dependence of the effective G-matrix. As shown by Jeukenne et al. [5|, in order to predict reliably the s.p. properties 
within the Brueckner theory, one has to go beyond the lowest-order BHF approximation by considering the high- 
order contributions in the hole-line expansion of the mass operator. The second-order term MJ is called the Pauli 
rearrangement term and it is induced by the medium dependence of the G-matrix via the Pauli operator in the BG 
equation |4l| . The Pauli rearrangement effect of the G-matrix describes the influence of the ground state two- hole 
correlations on the s.p. potential [H, [4^. The ground state NN correlations have been investigated extensively in 
literature [50l - l52j . The Pauli rearrangement have been shown to play its role mainly in the low momentum region 
around and below the Fermi surface where the ground state two-hole correlations are most effective, and it weakens 
the momentum dependence of the s.p. potential especially around the Fermi surface. The effect of ground state 
correlations not only is essential for getting a satisfactory agreement between the predicted depth of the microscopic 
BHF s.p. potential and the empirical value [&] and for restoring the Hugenholtz-Van Hove theorem which is destroyed 
seriously at the lowest BHF approximation |41| . but also plays a crucial role in generating a nucleon self-energy to 
describe realistically the s.p. strength distribution in nuclear matter and finite nuclei below the Fermi energy [10|. 
According to Refs. d, [4l|, the Pauli rearrangement contribution is calculated as follows from the G-matrix: 



MJ((.-,u.) = ^Mf'' = 



a'k'>kZ 



E 

ki<kZ 



E 

k2<kl 



|(fcfc'|G^^^'[e^(fci 



'(fc2)]|fclfe)P 



uj + e'-'ik') -e^(fci) -e^'(fc2) - ir] 



(5) 



Due to the NN correlations in asymmetric nuclear medium, the neutron and proton Fermi seas are partially depleted, 
and consequently the correlated momentum distributions differ from the uncorrelated ones. The third-order term 
in the hole-line expansion of mass operator is called the renormalization contribution and it takes into account the 
above effect of the depletion of the Fermi seas. The renormalization term MJ is given by [1, [2^, |4l| 



M3^(fc,a;) = ^^2 ih'){kh'\G^^' [uj + e^' {h'))\kh')A 



(6) 



where h' refers to the hole state with momentum smaller than fcp, and 



«2 {h') = - 



(7) 



is the depletion of the neutron or proton Fermi sea at the lowest-order approximation in asymmetric nuclear matter 
[E, EH, i.e., K2 {h') is the probability that a neutron or proton hole-state {\h'\ < kp) is empty. Hereafter, we shall 



use h and h' to denote the s.p. hole states below the Fermi momentum. As shown in Ref. [29j, it is an satisfactory 
approximation to replace in Eq. ([6]) the depletion coefficient ^2 (h') by its value at the averaged momentum inside 
the Fermi sea, i.e., = {h' = 0.75A;J.). We have AfJ(fc,a;) ~ — Ml'^ {k,uj). By taking into account the 



renormalization term MJ(fc,w), one may get the renormalized BHF approximation for the mass operator [41|, i.e.. 



Mlik, uj) = Ml{k, uj) + M^{k, uj) 



{k,Lu) 



(8) 



Similarly, one may consider a renormalization correction from the four hole-line terms to the second-order term Af J in 
order to take into account the fact that the hole-state k2 in Eq. ([5]) is partially empty ( see also Ref. ^ for symmetric 
nuclear matter ). Accordingly we obtain the renormalized M2, which is approximately given by (4]| : 



M^{k,uj) 



Mr 



{k,uj) 



(9) 



In terms of the ofF-energy-shell mass operator, one can readily calculate the neutron and proton momentum distribu- 
tions in asymmetric nuclear matter below and above the corresponding Fermi momentum [29| : 



dU[{k,uj) 



for k < kl 



(10) 



(jj=e^(k) 
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for fc > fcj, (11) 
where and L/J denote the real parts of and AfJ, respectively. 



duj 



III. RESULTS AND DISCUSSIONS 



In Fig. [T]we display the TBF effect on the predicted momentum distributions below and above the corresponding 
Fermi momenta in symmetric nuclear matter (/3 — 0) at two typical densities p = 0.17 and 0.34 fm~^, respectively. 
In the figure, the solid lines correspond to the results obtained by including the TBF; the dashed ones are calculated 
by adopting purely the AF18 two-body force alone. It is clear from Fig. [T]that, due to the many-body correlations 
induced by the NN interaction, the s.p. states below the Fermi surface are partly empty, and those above the Fermi 
surface can be partly occupied in the correlated ground state of nuclear matter. In the case of excluding the TBF, 
the density dependence of the momentum distribution, as a function of the ratio k/kp, is shown to be quite weak 
in the density region considered here, which is in good agreement with the previous EBHF calculation in Ref. [2^ 
by adopting the separable AVIA interaction and the prediction reported in Ref. [SS] by using the Green function 
method. One may notice that the TBF effect is negligibly small at low densities around and below the empirical 
saturation density po — 0.17 fm~^ of nuclear matter. This is consistent with the conclusion in Ref. {39| within the 
correlated basis function approach by adopting the Urbana wl4 interaction plus an effective three-body interaction. 
At the high density p — 0.34 fm~^ which is well above the saturation density, the TBF effect turns out to become 
noticeable. By comparing the solid curves and the corresponding dashed curves in Fig. [TJ it is seen that the TBF 
effect is to enhance the depletion of the momentum distribution below the Fermi momentum at high densities, i.e., to 
reduce the occupation probability of the hole states. This is readily understood since inclusion of the TBF is expected 
to induce stronger short-range correlations in dense nuclear medium as compared with the case of excluding the TBF. 
In both cases of including and excluding the TBF, the obtained depletions of the hole states at zero momentum 
(k = 0) are roughly 15% around the saturation density, compatible with the previous predictions reported in Refs. 
|26l [29i . [sgI [39I lioj . Since the depletion of the s.p. hole states well below the Fermi momentum in nuclear matter 
can be identified with the depletion of the occupation of the deeply bound s.p. levels in finite nuclei fs). [s^ . lioj . the 
present obtained depletion at saturation density is also consistent with the experimental result in Ref. |15| . 
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FIG. 1: (Color online) TBF Effect on the nucleon momentum distribution in symmetric nuclear matter (/3 = 0) for two densities 
0.17 fm"^ (left panel) and 0.34 fm"^ (right panel). 

The isospin T — tensor component of the N N interaction is expected to be crucial for determining the isospin 
dependence of the equation of state and s.p. properties of asymmetric nuclear matter 0, 0, Sll, [s^. The study of 
the neutron and proton momentum distributions as well as their isospin-asymmetry dependence may be helpful for 
understanding the properties of the short-range and tensor correlations in nuclear many-body systems 2^ 36]. In 
Fig. [2]we report the neutron and proton momentum distributions below and above their respective Fermi momenta in 
asymmetric nuclear matter at various asymmetries /3 = 0, 0.2, 0.4, 0.6, and 0.8 for two typical densities p = 0.17 and 
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0.34 fm"'^, respectively. In the figure, the results are obtained by adopting purely the two-body interaction and 

the TBF is not included. It is clearly seen that the neutron and proton momentum distributions in asymmetric nuclear 
matter are different from their common distribution in symmetric nuclear matter and the depletions of the neutron 
and proton hole states depend sensitively on the isospin-asymmetry. As the isospin-asymmetry j3 increases, the 
occupation probability of the neutron hole states below the neutron Fermi sea becomes larger while the occupation 
of the proton hole states gets smaller with respect to their common values in symmetric nuclear matter; that is, 
increasing the asymmetry leads to a reduction of the depletion of the neutron hole states, while it enhances the 
depletion of the proton hole states. The above result has also been found in Ref. [s^ within the framework of the 
Green function approach. Such an isospin-asymmetry dependence of the neutron and proton momentum distributions 
in asymmetric nuclear matter implies that at a higher asymmetry, the effect of the short-range correlations, induced 
by the TV TV interaction, becomes stronger (weaker) on protons (neutrons), and can be understood according to the 
isospin-asymmetry dependence of the effect of the tensor component of the TV TV interaction in asymmetric nuclear 
matter. As is well known, the isospin T = SD tensor component may induce a strong short-range correlation in 
nuclear medium and it plays a dominant role in determining the isospin vector parts of the properties of asymmetric 
nuclear matter fZl], Issl - lssj . As the neutron excess increases, the effect of the T = SD tensor channel on protons 
(neutrons) from the surrounding neutrons (protons) becomes stronger (weaker). Accordingly the SD tensor component 
is expected to induce a larger (smaller) depletion of the proton (neutron) Fermi sea at a higher asymmetry. The definite 
evidence for the strong enhancement of the np short-range correlations over the pp and nn correlations observed at 
JLab 22J has provided an experimental indication for the dominant role played by the short-range tensor components 
of the TV TV interactions in generating the TV TV correlations [1^ . The importance of the tensor force in determining the 
isospin-asymmetry dependence of the neutron and proton momentum distributions in asymmetric nuclear matter has 
also been confirmed in Ref. [36| where the authors have calculated the depletion of the nuclear Fermi sea by adopting 
various TV TV interactions and they find that the iso-depletion [defined as the difference of the neutron and proton 
occupation of the lowest momentum state in asymmetric nuclear matter, i.e., n"(A; = Q) — •nP{k = 0)] obtained by 
the AVA' potential which has no tensor component is much lower than the predictions by the other TV TV interactions 
including the tensor components, especially at suprasaturation densities. 
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FIG. 2: (Color online) Neutron and proton momentum distributions in asymmetric nuclear matter at various asymmetries for 
two densities 0.17 fm~'^ (left panel) and 0.34 fm^'^ (right panel). The results are calculated without including the TBF. 

In Fig. [21 we show the neutron and proton momentum distributions below and above their respective Fermi mo- 
menta, predicted by adopting the ^VIS interaction plus the TBF, in asymmetric nuclear matter at various asym- 
metries /3 = 0, 0.2, 0.4, 0.6, and 0.8 for two typical densities p = 0.17 and 0.34 fm~'^, respectively. It is seen that, 
after including the TBF in the calculation, the predicted isospin-asymmetry dependence of the neutron and pro- 
ton momentum distributions remain the same as that obtained by adopting purely the ^^18 two-body force. By 
comparing the corresponding curves in Fig. [5] and Fig. [31 one may notice that the TBF effect is negligibly small at 
sub-saturation densities. At densities well above the saturation density, the TBF leads to an overall enhancement 
of the depletions of both the neutron and proton distributions below their respective Fermi momenta since it may 
generate extra short-range TV TV correlations which become strong enough at high densities. 

In order to see more clearly the isospin dependence and the TBF effect, in Fig. [31 we display the proton and neutron 
momentum distributions at zero momentum fc = as functions of the isospin-asymmetry 13 in the two cases with (solid 
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FIG. 3: (Color online) The same as Fig. [2l but the results are obtained by including the TBF. 
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FIG. 4: (Color online) Neutron and proton momentum distributions at zero momentum in asymmetric nuclear matter vs. 
isospin asymmetry j3 for two densities 0.17 fm~"^ (left panel) and 0.34 fm~"^ (right panel). The results are obtained for the two 
cases of including the TBF (solid curves) and excluding the TBF (dashed curves). 



curves) and without (dashed curves) including the TBF. It is clearly seen that the proton (neutron) occupation of the 
zero momentum state decreases (increases) almost linearly as a function of asymmetry, indicating the short-range and 
tensor correlations become stronger (weaker) for protons (neutrons) at a higher asymmetry in neutron-rich nuclear 
matter. A similar result has been reported in Ref. [13] for the averaged neutron and proton occupation probabilities 
below their respective Fermi seas at the empirical saturation density. The above predicted isospin splitting [i.e., 
n"(A; = 0) > n^(fc = 0)] of the neutron and proton occupation of the lowest states in neutron-rich nuclear matter has 
also been found in Ref. [s^ by using the Green function approach. In Ref. [s^ , it is shown that at finite temperature 
the thermal effect may destroy the linear dependence of n"(fc = 0) and nP{k = 0) on asymmetry /3 at high enough 
asymmetries. As expected, around and below the saturation densities, the TBF effect on the proton and neutron 
momentum distributions is negligibly small in agreement with the result of Ref. "sO*]. However, at high densities ( 
for example p — 0.34 fm~'^ ) well above the saturation density, the TBF may affect sizably the neutron and proton 
occupations of their lowest states. Inclusion of the TBF leads to an overall enhancement of the depletion of the 
neutron and proton hole states. It is noticed from Fig. 2] that the TBF has almost no effect on the iso-dcpletion 
n"(0) — nP{0) up to /o = 0.34fm^'^, which provides a support for the result of Ref. [3^ where it has been shown that 
once the tensor components fitting the experimental phase shifts are included, various modern NN interactions lead 
to almost the same iso-depletion although the momentum distributions of neutrons and protons predicted by those 
different A'^A'^ interactions can be quite different. 

Before our summary, we shall give a brief discussion about the possible relevance of the present results for the 
dense and highly asymmetric nuclear matter in the interior of neutron stars, especially the cooling of neutron stars 
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via neutrino emission and the nucleon pairing in neutron stars. Let us consider a neutron star consisting of neutrons, 
protons and electrons at /3-equilibriuni, i.e., a (n,p, e) neutron star model. For a {n,p,e) neutron star, there are two 
different kinds of URCA processes for neutrino emission. One is the direct URCA process, the other is the modified 
URCA process. The direct URCA process may lead to a much faster cooling of neutron stars than the modified 
URCA process. Under the assumption of a free Fermi momentum distribution of nucleon (i.e., the Fermi sea is fully 
occupied and the particle states above the Fermi momentum are completely empty) , the direct URCA process is only 
allowed if the proton/neutron ratio x = Pp/pn is greater than a threshold value of Xc = 1/8 [which corresponds to a 
proton fraction of Yp = Pnjipn + Pp) = 1/9 and an isospin-asymmetry of /3 = 7/9 ~ 0.8] in order to guarantee the 
momentum conservation [56| . The strong depletion of the proton Fermi sea and the partial occupation of the proton 
states well above the Fermi momentum, induced by the short-range np correlations in dense and highly asymmetric 
nuclear matter, may affect considerably the direct URCA process as has been discussed in detail by Frankfurt et al. in 
Ref. [24|| where it is shown that the modification of the proton momentum distribution in neutron star matter due to 
the short-range correlations leads to a significant enhancement of the neutrino luminosity of the direct URCA process 
for temperatures much less than 1 MeV, and the direct URCA process may even have probability to occur for a small 
value of X < 0.1. Nucleon pairing in asymmetric nuclear matter plays an important role in determining the cooling 
rate of neutron stars (57| and its strength has been shown to be rather sensitive to the medium effect induced by 
nucleon- nucleon correlations (ssj . The large depletion of the proton Fermi sea in highly asymmetric nuclear matter due 
to the T = np short-range correlations, obtained in the present calculation, tends to reduce considerably the proton 
pairing in neutron stars [24| . Another interesting topic is the thermal transport parameters in dense asymmetric 
nuclear matter which have special importance for the damping of nonradial modes of neutron stars and may depend 
sensitively on the depletion of nucleon distribution [59l - [6^ . In general, the transport parameters, including the shear 
viscosity and thermal conductivity, may be calculated within the framework of the Landau Fermi liquid approach 
in which a free Fermi momentum distribution is assumed as the equilibrium distribution [59| . Therefore, a strong 
depleted proton Fermi sea in dense and supra-dense asymmetric nuclear matter is expected to affect substantially the 
transport properties in the interior of neutron stars. 

IV. SUMMARY 

In summary, we have investigated the TBF effect on the proton and neutron momentum distributions in asymmetric 
nuclear matter within the framework of the EBHF approach by adopting the AU18 two-body interaction supplemented 
with a microscopic TBF. In symmetric nuclear matter {(3 = 0), the obtained depletion of the hole states deep inside 
the Fermi sea is roughly 15% at the empirical saturation density, in agreement with the previous predictions [26. . ,29l . 
m, [1^ and the experimental value in Ref. [l^. In asymmetric nuclear matter {(3 > 0), the neutron and proton 
momentum distributions turn out to become different and may split with respect to their common distribution in 
symmetric nuclear matter. It is shown that increasing the isospin-asymmetry /3 tends to enhance the depletion of the 
proton Fermi sea while it reduces the depletion of the neutron Fermi sea, which implies that at a higher asymmetry, the 
effect of the tensor correlations induced by the NN interaction may become stronger on protons while it gets weaker 
on neutrons. At zero momentum, the neutron occupation probability increases while the proton occupation decreases 
almost linearly as a function of asymmetry. The present obtained isospin dependence of the neutron and proton 
momentum distributions in asymmetric nuclear matter is in good agreement with the recent prediction in Ref. [36| 
within the framework of the Green function method. At low densities around and below the nuclear saturation 
density, the TBF effect on the predicted momentum distributions is found to be negligibly weak in agreement with 
the conclusion of Ref. [39[ . At high densities well above the saturation density, the TBF is expected to induce strong 
enough extra short-range correlations and its effect turns out to become noticeable. In dense asymmetric nuclear 
matter, inclusion of the TBF effect may lead to an overall enhancement of both the depletion of the neutron and 
proton Fermi seas for all the asymmetries considered. Although the TBF affects sizably the neutron and proton 
momentum distributions at high densities well above the saturation density, its effect on the iso-depletion of the 
nuclear Fermi sea (i.e., the difference of the neutron and proton occupation probabilities) in asymmetric nuclear 
matter is shown to be quite small in the density region up to two times saturation density. The present results are 
expect to have significant implication for the cooling, nucleon pairing and transport properties of neutron stars. 
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